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Abstract: Using the Resonance Chiral Theory Lagrangian, we perform a calcu- 
lation of the vector form factor of the pion at the next-to-leading order (NLO) 
in the l/Nc expansion. Imposing the correct QCD short- distance constraints, one 
fixes the amplitude in terms of the pion decay constant F and resonance masses. 
Its low momentum expansion determines then the corresponding 0(p 4 ) and 0(p 6 ) 
low-energy chiral couplings at NLO, keeping control of their renormalization scale 
dependence. At /io = 0.77 GeV, we obtain L 9 (/i ) = (7.9 ± 0.4) • 10~ 3 and 
C 88 (/io) - C 90 (/i ) = (-4.6 ± 0.4) • 10" 5 . 
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1. Introduction 

Effective field theories (EFT) are nowadays the standard tool to investigate the low- 
energy dynamics of Quantum Chromodynamics (QCD). In particular, the chiral sym- 
metry is a crucial ingredient for the understanding of the light quark interactions. The 
dynamics of the pseudo-Goldstone bosons from the spontaneous symmetry breaking 
is provided by the corresponding EFT, Chiral Perturbation Theory (xPT), with a 
perturbative expansion in powers of light quark masses and external momenta [jl], . 
This allows a systematic description of the long-distance regime of QCD, at energies 
below the lightest resonance mass. The precision required in present phenomeno- 
logical applications makes necessary to include corrections of 0(p 6 ). While many 
two-loop ;\FT calculations have been already carried out ||, the large number of 
unknown low-energy constants (LECs) appearing at this order puts a clear limit to 
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the achievable accuracy. The determination of these couplings is compulsory to 
achieve further progress in our understanding of strong interactions at low energies. 

In the resonance region, E ~ Mr, the chiral counting breaks down and the new 
heavier degrees of freedom -the resonances- have to be explicitly incorporated into 
the theory. A suitable alternative is then provided by the 1/Nc expansion in the limit 
of a large number of colours, N c — >• oo [|, |], ||. Assuming confinement, the strong 
dynamics is given at large Nq by tree-level diagrams with an infinite number of pos- 
sible hadronic exchanges. This corresponds to the tree approximation of some local 
Lagrangian, being meson loops suppressed by higher powers of l/N c j4|. Resonance 
Chiral Theory (R%T) provides an appropriate framework to incorporate these mas- 
sive mesonic states within a chiral invariant phenomenological Lagrangian [J7L H, 01. 
The operators of the R%T action are constructed such that they remain unchanged 
under flavour transformations U(3)l ®U(3)r. After integrating out the heavy fields, 
the xPT Lagrangian is recovered at low energies with explicit values of the chiral 
LECs in terms of resonance parameters. The short- distance properties of QCD impose 
stringent constraints on the R;\T couplings and provide important information for 
the extraction of the low-energy %PT parameters. The amplitudes are thus enforced 
to follow the known high-energy QCD behaviour, introducing in the long-distance 
description important information from the underlying theory [[5], 

Clearly, we cannot determine at present the infinite number of meson couplings 
which characterize the large-iVc Lagrangian. However, one can perform useful ap- 
proximations in terms of a finite number of meson fields. Truncating the infinite tower 
of mesons to the lowest resonances with h , ++ , 1 and 1 ++ quantum numbers, 
one gets a very successful prediction for the 0{p A ) couplings at large Nc |J. 
Already at this level the comparison with experimental determinations of the 0(p 4 ) 
chiral couplings shows a remarkable agreement. Some 0(p 6 ) LECs have been also 
estimated in this way, by studying appropriate sets of Green functions (see ref. PJ 
and references therein). All the required terms in the R%T Lagrangian that may 
contribute to the 0(p 6 ) LECs at LO in 1/N C were classified in ref. PI. 

Since chiral loop corrections are of next-to-leading order (NLO) in the 1/Nc 
expansion, the large-iVc determination of the LECs is unable to control their 
renormalization-scale dependence. First analyses of resonance loop contributions to 
the running of L w (fi) and L 9 (/x) were attempted in refs. [|K| and []TT[] , respectively. 
In spite of all the complexity associated with the still not so well understood renor- 
malization of R%T [JTT], [12], [13|, [14], |5|, [L6[], these pioneering calculations showed the 
potential predictability at the NLO in 1/Nc- 

Using dispersion relations we can avoid the technicalities associated with the 
renormalization procedure jTS], [T7|, [TJ|. This allows one to understand the underly- 
ing physics in a much more transparent way. Still, a fully equivalent diagrammatic 
calculation is possible, although the derivation and presentation is slightly more cum- 
bersome |10|, O, pll. In particular, the subtle cancellations among many unknown 
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renormalized couplings found in ref. [IT] and the relative simplicity of the final result 
can be better understood in terms of the imposed short- distance constraints within 
the dispersive approach. Following these ideas we determined, up to NLO in 1/Nq, 
the couplings L 8 (fi) and C^fi) in ref. |17| and Lio(fi) and C%i(n) in ref. ff8|. In this 
article we present the study of the vector form factor (VFF) of the pion, which allows 
us to estimate the coupling L 9 (/x) and the 0{p & ) combination C 88 (/i) — Cgo(^) 
up to NLO in l/N c . 

In order to establish the notation, the R%T Lagrangian is introduced in the next 
section. The analysis of the VFF in the resonance region is performed in section |3], 
while section f| contains the determination of L 9 (fi) and C 88 (/i) — C 90 (/i). A summary 
of our results is finally given in section |5|. In order to ease the reading of the text, we 
have shifted the technical details on the calculation of the spectral function, the full 
VFF and the chiral coupling expressions to the Appendices. 



2. The Lagrangian 

We will adopt the Single Resonance Approximation (SRA), where just the lightest 
resonances with non-exotic quantum numbers are considered. 1 On account of the 
large- N c limit, the mesons are put together into U(3) multiplets. Hence, our degrees 
of freedom are the pseudo-Goldstone bosons (the lightest pseudoscalar mesons) along 
with massive multiplets of the type V(l~~), A(l ++ ), S(0 ++ ) and P((T + ). With them, 
we construct the most general action that preserves chiral symmetry. Since we are 
interested in determining the ;\RT low-energy constants and the study of the short- 
distance behaviour, the chiral limit will be taken all along the paper. No information 
is lost as the chiral LECs are independent of the light quark masses. 

Resonance Chiral Theory must satisfy the high-energy behaviour dictated by 
QCD. To comply with this requirement we will only consider operators constructed 
with chiral tensors of 0(p 2 ); interactions with higher-order chiral tensors tend to vi- 
olate the asymptotic short- distance behaviour prescribed by QCD || Likewise, 



it has been shown in some cases that resonance operators with higher number of 
derivatives can be simplified into terms with less derivatives, terms without reso- 
nances and operators that contribute to other hadronic amplitudes, by means of the 
equations of motion and convenient meson field redefinitions 0, |9|, [TT], [12], 13], [19 . 



1 In ref. PQ| , it has been argued that large discrepancies may occur between the values of the 
masses and couplings of the full large- Nc theory and those from descriptions with a finite number 
of resonances. Even in this case, it is found that one can obtain safe determinations of the LECs as 
far as one is able to construct a good interpolator that reproduces the right asymptotic behaviour 
at low and high energies. Further issues related to the truncation of the spectrum to a finite number 
of resonances are discussed in ref. lEli. 
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The different terms in the Lagrangian can be classified by their number of reso- 
nance fields: 



Ri R±,R2 



where the dots denote operators with three or more resonance fields, and the indices 
Ri run over all different resonance multiplets, V, A, S and P. The term with only 
pseudo-Goldstone bosons is given by [H 



£c = j(t»/ + X+)- (2-2) 



The second term in eq. ( |2.1|) corresponds to the operators with one massive reso- 
nance JF|, 

Cv = ^ {V - R '" > + ''iM {v ^ u "' u " ]) ' 

C s = c d (Su^ 1 ) + c m {Sx+) , 

C P = id m (Px-) ■ (2.3) 

The Lagrangian contains the kinetic resonance terms and the remaining oper- 

ators with two resonance fields [0, [TTJ] . We show only the terms that contribute to 
the vector form factor of the pion, taking into account that here we just consider the 
lowest-mass two-particle absorptive channels, with two pseudo-Goldstone bosons or 
one pseudo-Goldstone and one resonance. In the energy range we are interested in, 
exchanges of two heavy resonances are kinematically suppressed. Hence, the relevant 



operators are 




AC S A = 


Af A ({V^,^K), 


AC S p = 


Af> a {V a S,P}>, 


AC PV = 


zaH[v^p,^k>, 


ACva = 


i\l A ([V^,A va \hl) + 



>VA, 



+ i\Z A ([V a V^,A™]u») +i\v A ([V a V lxu ,A^]u a ) . (2.4) 

All coupling constants are real, the brackets (...) denote a trace of the corresponding 
flavour matrices, and the standard definitions for the u^, x±-> f± an d h^ v chiral 
tensors of pseudo-Goldstones are provided in refs. f7|, || . 
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Our Lagrangian C-r x t satisfies the N c counting rules for a theory with £7(3) 
multiplets. Therefore, only operators that have one trace in the flavour space are 
considered. Note that local terms with two traces in flavour space, which are of 
NLO in 1/Nc, cannot contribute at tree-level to the VFF because the final two-pion 
state has isospin 1 = 1. The different fields, masses and momenta are of 0(N C ) in 
the 1/Nc expansion. Taking into account the interaction terms, one can check that 
F, Fy, Gy, Fa, Cd, c m and d m are 0(y/Nc) and the Af li?2 are 0(N C ). The mass 
dimension of these parameters is [F] = [Fy] = [Gy] = [Fa] = [q] = [c m ] = [d m ] = E 
and \'< "■' /•;". 

Note that the U(3) equations of motion have been used in order to reduce the 
number of operators. For instance, terms like (PV^u^) are not present in eq. (|2.3|), 



since they can be transformed into operators that, either have been already consid- 
ered, or contain a higher number of mesons by means of the equations of motion and 
convenient meson field redefinitions 0. 

The R%T Lagrangian fl2.1|) contains a large number of unknown coupling con- 
stants. However, as we will see in the next section, the short-distance QCD constraints 
allow us to determine many of them. In the observable at hand and with our assump- 
tions, we initially have ten couplings or combinations of them (F, Fy, Gy, Fa, q, 
Xf A , Xf p , Af y , -2A^ A + A^ A and 2A^ A - 2A^ A + Aj A + 2A^ A ) and four resonance 
masses (My, Ma, M s and M P ). As we will see in section |3|, after imposing a good 
short-distance behaviour of this observable, the number of parameters reduces to 
three couplings (F, Gy and Fa) and three masses (My, Ma and Ms)- The Weinberg 
sum-rules associated with the left-right correlator allow us to further reduce the 
number of inputs; the amplitude is finally determined in terms of just F and the three 
masses My, Ma and Ms- The role of the information coming from the underlying 
theory is thus fundamental. 

3. The vector form factor of the pion 

Our observable is defined through the two pseudo-Goldstone matrix element of the 
vector current: 

( * + ( Pl ) 7r-(p 2 ) | i (u^u - dYd) |0 ) = F(s) ^ - p 2 r , (3.1) 

where s = (p± +P2) 2 - At very low energies, F(s) has been studied within the 
framework up to 0(p 6 ) [0, |23|]. R%T and the l/N c expansion have also been used 
to determine F(s) at the p meson peak, including appropriate resummations of sub- 
leading logarithms from two pseudo-Goldstone channels [^4], |25|. A first systematic 



study of the VFF at NLO in 1/Nc was performed in ref. |TT|]. Although the gen- 
eral structure was well established there, the present article answers and solves three 
important questions raised in that previous paper: 
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Figure 1: Tree-level contributions to the vector form factor of the pion. A single line 
stands for a pseudo-Goldstone boson while a double line indicates a resonance. 



In ref. JTTJ] only operators with at most one resonance field were included (except 
for the kinetic resonance terms) []7J]. However, as suggested in the appendix C 
of that article, this assumption is not really justified and leads to problems 
with the asymptotic short distance behaviour. In the present paper, we have 
considered all the operators needed to describe the absorptive cuts with two 
chiral pseudo-Goldstones and those with one pseudo-Goldstone and one reso- 
nance, being higher thresholds with two resonances highly suppressed in the 



energy region that we consider [[18 



Due to this first issue, in ref. |TT] the logarithmic part of J~(s) was badly 
behaved at high energies. It was not possible to enforce a vanishing form factor 
at s — > oo without the inclusion of new hadronic operators in the leading 
Lagrangian. The inclusion of those terms in the present article will allow us to 
recover the expected high-energy dependence for the VFF in QCD [|26| . 



The final result of ref. [|11|] contained the unknown R/\T couplings L$ and 
C$s — Cgo, which are the analogous ones to the xPT LECs L 9 and C$s — Cgo- 
In the present work, they are fully determined by means of the high-energy 
matching with QCD [Tb 



Within Resonance Chiral Theory the diagrams contributing to the VFF at lead- 
ing order in 1/Nc are shown in figure |l[ They generate the result 

= i + ^5^7- (3-2) 

Considering that the form factor is constrained to be zero at infinite momentum 
transfer [E5|, the vector couplings should satisfy 

FyGy = F 2 , (3.3) 



which implies 



M 2 
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Figure 2: One-loop contributions to the vector form factor of the pion with absorptive 
cut. A single line stands for a pseudo-Goldstone boson while a double line indicates a 
resonance. 



The subleading corrections can be calculated by means of dispersive relations. 
Once the one-loop absorptive parts of Tr x t (s) are known, one can reconstruct the 
full form factor up to appropriate subtraction terms. We can separate then the leading 
and subleading parts of the amplitude in the form 

M 2 

FrxT (s) = M 2 V _ S + jr ( s ) NL o , ( 3 - 5 ) 

with .F(s) NLO containing the one-loop contribution and the subleading part <5 NLO of 
the resonance coupling combination FyGy j ' F 2 = 1 + 5 NLO (for details see appendix 

@): 

^(s) NLO = -Tfj— + ^(s) . (3.6) 

My — S 

The explicit form for the subtracted one-loop amplitude JF {s) can be found in 
Appendices [A| and 0, being fully determined by the spectral function ImJ-'(s) through 
a once-subtracted dispersion relation. It vanishes at s = and has no contribution to 
the real part of the pole at s = My. The subleading correction to the couplings, £ NLO , 
is fixed by means of the high-energy matching after demanding that it cancels the 
bad behaviour of J zU {s) = 5 NLO + 0{s~ x ) when s — > oo. Furthermore, the NLO term 
J-"(s) NLO can be neatly separated into its different contributions from the various two- 



NLO 

,,. LO can be neatly separ; 
meson absorptive channels J 7 (s) NLO |mi,m 2 5 given by the corresponding F u {S)\ mi , m2 
and the consequent 8 NLO \ mi)m2 . These details are relegated to Appendices and 0. 
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Although in this article we follow the procedure of refs. JIT], [RJ, our results can be 
also derived in an utterly equivalent way through a Feynman diagram computation 
and the standard renormalization procedure. This derivation is slightly more complex 
and its detailed explaination is relegated to appendix ||. 

We will consider only the effects of absorptive loops with two pseudo-Goldstones 
(nn) or with one pseudo-Goldstone and a resonance (Rtt). Two-resonance channels 
RR' have their thresholds at (Mr + M' R ) 2 > 2 GeV 2 and their impact on the LEC 



determination is expected to be negligible [fT8[- Taking this into account, we extract 



our R/\T form factor through the following short- distance matching procedure: 

1. Determine the spectral function of the considered absorptive cuts (tttt and Rir). 
The full expressions are shown in eqs. ( |B.1|) , (|B.2|) and (|B.3|) of appendix |B[ 



2. We demand ImJ-'(s) to be well-behaved at high energies, i.e., it must vanish 
when s — > oo. In the present work, we will actually impose this constraint 
channel by channel, i.e., we will demand that each separate two- meson cut 
ImJ^s)!^ ^ vanishes at s — > oo. For spin-0 mesons this must be so as its one- 
loop contribution to the spectral function is essentially its VFF at LO (which 
vanishes at infinite momentum) times the partial-wave scattering amplitude at 
LO (which is upper bounded). For higher spin resonances the derivation is more 
cumbersome as the Lorentz structure allows for the proliferation of form factors 
and the unitarity relations are not that simple. Still, in many situations it has 
been already found that amplitudes with massive spin-1 mesons as final states 
must go to zero at high energies even faster, due to the presence of extra powers 
of momenta in the unitarity relations coming from intermediate longitudinal 
polarizations [0. In summary, we will assume ImF^s)]^^ — > when s — > oo 
for every absorptive two-meson cut under consideration, regardless of the spin 
of the intermediate mesons. 

In the case of the nn cut we have found two constraints, which are consistent 
with the literature, 

F V G V = F 2 , 3G 2 v + 2c 2 d = F 2 , (3.7) 

where the first one coincides with eq. ( |3.3|) , that is, with the constraint obtained 
with the vector form factor at leading-order ||. The second one was derived 



in ref. from the LO im scattering amplitude. It is interesting to remark 
that the ca = limit of this second relation, Gy = Fj a/3, has been obtained 
recently from a study of t~ — > P~ 7 v T decays [P = n, K) p5| . We have used 
these constraints to fix Fy and c 2 d . 

For the Pn cut, the only possible solution is to kill the whole contribution by 
means of 

Af v = , (3.8) 
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which is consistent with the large- A*^ constraint from the vector form factor 
into Pn, studied in ref. 



The analysis of the An cut leads to more than one real solution. We have 



chosen the solutions consistent with previous works [|15|, |I8|, where the NLO 
contributions in 1/Nc to the Uyv{s) correlator coming from tree-level form 
factors to resonance fields were studied: 

-2A^ A + Af = , - Af + \J A + 2Af = ^ , 

ASA = F A G V {M\ — 4 My) 
1 3V2Mlc d F v 

The first two constraints, in the first line, come from the analysis of the 
An vector form-factor. The last relation with Af A is then needed to make 
ImJ^s)!^ — )• for s — )• oo. 

After imposing the relations (^.7[), ( |3.8| ) and ( PT5| ) the spectral functions can 
be expressed in terms of Gy, Fa, F and masses, as shown in eqs. ( |B.6[ ), ( |B.7j ) 
and (gj ). 

3. The spectral function is now ready for the once-subtracted dispersion relation 
provided in the appendix |A] in eq. ( |A.4| ) , which allows to reconstruct the full 



form factor up to the pole position at s = M v and the real part of its residue. 

Finally, we impose that the whole FrxtIs) vanishes at short distances -not 
only its imaginary part-. This fixes the real part of the residue at s = M, 2 



v 



and, consequently, the NLO correction <5nlo i n eo ±- flj-1)- m order to ease the 
reading of the manuscript, the complicated expressions for the well-behaved 
contributions to the different channels are provided in appendix [C], in eqs. ( p.l|) , 
O and Q. 



4. The chiral couplings Lg(fi) and Cg$(fi) — Cgo(fi) 



The low- momentum expansion of J-'(s) is determined by %PT ||, The corre- 
sponding expression in the chiral limit reads 



X PT 



1 + 



2s 

V 2 



LM + 



32tt 2 



log- 



4s 2 



C 88 (/i) - C 90 (/i) 



r (L) _ r (L) 
1 88 1 90 

327T 2 



5 -s 
3- l0g ^ 



with 



(4.1) 



1 

4' 



r (i) _ r 

1 88 1 90 



(L) 



2Lx 



L2 
3 



2 4 



9 



(4.2) 
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The couplings F 2 , Lg, C 88 /F 2 and C 90 /F 2 are of 0(N C ), while r 9 , T^/F 2 and 
TfJ/F 2 are of 0{N C ) and represent a NLO effect. 

The low-energy expansion of eqs. ( |3.4| ) and (|3.5p, obtained, respectively, within 
Resonance Chiral Theory at leading-order and at next-to- leading order in the 1 /Nc 
expansion, allows to determine the chiral couplings Lg and C 88 — C90 & t LO and at 
NLO. 

4.1 The large- N c limit 

At leading-order in 1/Nc, eq. (|L1|) becomes 



J- x pt (a) = 1 + {^9 + O(N c ) } - ^ {C 88 - C 90 + 0(N C )} + 0{s 3 ) . (4.3) 
Within RxT in the large- Nc limit, eq. ( |3.4[) can be now expanded at low energies: 

The matching between (fO|) and (fL4|) fixes L 9 and C 88 - C 90 at LO g 

F 2 F 4 

4.2 L 9 (//) and C 88 (//) - C 90 (//) at NLO 

Following the same steps as before, let us determine the related 0(p 4 ) and 0(p 6 ) low- 
energy constants by matching eq. ( f4.1|) and the low-energy expansion of eq. fl3.5|), 

r 9 



f, 



RxT 




where the £( 2n ) are the relevant 0(s n ) coefficients of the low-energy expansion of 
J NLO (s), once the structure coming from the one- loop diagram has been sub- 
tracted from the tttt channel. The separated contributions fm"m. 2 from each absorp- 
tive two-meson cut F NLO (s)\ mi , m2 are provided in appendix [D|, being each of them 
independent of the renormalization scale \i. 

By comparing the xPT expression (|4.1| ) to the R%T low-energy expansion (|4.6|) , 
it is straightforward to estimate the chiral LECs Lg(fi) and C 88 (/x) — C 90 (fi) up to 
NLO in l/N c : 

T , % f 2 F<2) r 9 My 

F 4 r (L)_ r (L) ,, 2 l 4 -0 



-10- 



where 



rW _ p W 3Gy _ c 2 , -F\/GV F 2 - 3Gy _ F 2 + 3Gy 

ss 90 ~ 8M 2 4M 2 + 8M 2 ~ 8M| 8M 2 1 



matches the corresponding 0(p 6 ) running at NLO in 1/Nq. Note that the large- Nc 

CP" 3G^ TP Cr 

relations L 2 = 2L\ = ^r, L 3 = —j^r + jm? an d L 9 = |7| have been used in 

eq. (|4.2|) . The high-energy constraints F V G V = F 2 and 2c 2 , = F 2 — 3Gy of eq. flvT 



have been employed to obtain the result on the r.h.s. of eq. (|4.8j ). 
4.3 Phenomenology 

Using My ~ 0.77 GeV and F ~ 89MeV, one gets the large- Nc estimates from 
eq. (|Op: L 9 ~ 6.7 • 1(T 3 and C 88 - C 90 ~ -4.5 • 10~ 5 . At /i = 770 MeV, the 
phenomenological determinations L 9 (/i ) = (6.9 ± 0.7) ■ 10" 3 [@, || and L 9 (fi ) = 
(5.93 ± 0.43) •10- 3 ,C r 88 (/io)-C9o(^o) = (-5.5 ± 0.5) • 10~ 5 ||, obtained respectively 
from an 0(p 4 ) and an (9(p 6 ) ChPT fit, agree approximately with the LO estimates. 

Large- Nc estimates are naively expected to approximate well the couplings at 
scales of the order of the relevant dynamics involved (// ~ Mr). However, they always 
carry an implicit error because of the uncertainty on //. This theoretical uncertainty 
is rather important in couplings generated through scalar meson exchange, such as 
L 8 (fi). In the present case, it also has a moderate importance. The size of the NLO 
corrections in 1/Nq to L 9 (fi) and C 88 (fi) — C 90 (/i) can be estimated by regarding 
their variations with fi. These are respectively given by 

9LM _ r 9 _ Q g io _ 3; a (CM - CM) _ rg } - iff ^ Q Q 1Q _ 5 ^ 

<91og/i 2 327T 2 d log fx 2 32tt 2 

(4.9) 

So far, we have been working within a U(3)l ® U(3)r framework, but we are 
actually interested on the couplings of the standard SU(3)l <8> SU(3)r chiral theory. 
Thus, a matching between the two versions of xPT must be performed. Nonethe- 
less, on the contrary to what happens with other matrix elements (e.g. the S — P 
correlator [17]]), the spin-1 two-point functions do not gain contributions from the 



C/(3)-singlet chiral pseudo-Goldstone; the rj\ does neither enter at tree-level nor in 
the one-loop correlators. Therefore, the corresponding LECs are identical in both 
theories at leading and next-to- leading order in 1/Nq: L 9 (/j,) u ^ = L 9 (fi) su ^ 3 \ 

(C w (jm) - c w Qj.)) m = (CmOO - CM) SU(3) . 

The needed input parameters are defined in the chiral limit. We take the 
ranges g ^ M v = (770 ± 5) MeV, M s = (1090 ± 110) MeV and F = (89 ± 2) MeV. 
The resonance couplings Gy and Fa can be fixed in terms of F and masses if one 
considers the short- distance conditions obeyed by the left-right correlator ||. The 
constraint of eq. ( |3.3|) , coming from the vector form factor of the pion, and those from 



the first and second Weinberg sum rules [£2| determine the vector and axial-vector 
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1st Approach 


2nd Approach 


10 3 • L 9 at LO 


6.68 


6.68 


io 3 • m 
10 3 • m 


0.11 
0.00 
1.12 


-0.04 
0.00 
1.00 


105 . (Co R - Con) at LO 


-4.46 


-4.46 


io 5 • 

io 5 • e2 


0.76 
0.00 
-0.88 


0.71 
0.00 
-0.73 



Table 1: Different contributions to the chiral couplings within the two numerical ap- 
proaches explained in the text. 



couplings at LO in l/N c [|T|, IS 



F 2 _ F 2 Ml _ 2 M\-M V 2 _ 2 My 

7 1 1 M A -Mf Gv ~ F M\ ' Fa ~ F M'i-Mf (410) 

with Ma > My. Due to the large width of the ai(1260) meson, the determina- 
tion of the Lagrangian parameter Ma is far from trivial. From the observed rates 
r(p° ->• e+e-) = (7.02 ± 0.13) keV |2p] and T ( fll -)> 7r 7 ) = (650 ± 250) keV and 
considering QODQ , one finds M A = (938 ± 13) MeV and M A = (960 ± 80) MeV. 



Another large-iVc 1 determination of Ma was obtained in ref. |30[ from the study of 
the 7i — > ez/ e 7 decay, which yields Ma = (998 ± 49) MeV. We cannot use the in- 
formation coming from T(p — > 2tt) = (149.4 ± 1.0) MeV [^9] in order to determine 
Ma, since Gy is constrained by eq. (|3.7|) to be smaller than F/v3, which results in 
Ma < 940 MeV. In spite of the dispersion of values for Ma, one gets a consistent 
description in the range Ma = (920 ± 20) MeV, which we will take as our input. The 
resulting numerical predictions for the LECs are 

Lg(//o) = (7.9±0.4)-10- 3 , 
C 8 sM - C 90 (/i ) = (-4.6 ± 0.4) • IO" 5 , (4.11) 



being /z the usual renormalization scale, yU = 770 MeV. 

Alternatively, one could also use the phenomenological values for Gy, Fa and 
the axial-vector mass, instead of fixing them through the Weinberg sum-rules. Thus, 
one may employ M A = (1200 ± 200) MeV |2|, and F A = (120 ± 20) MeV, from 
the observed rate T (ai — > rcy) = (650 ± 250) keV p9 |. The constraint of eq. (|3.7| ) 
implies that Gy < Fj v^, so that we take the range Gy £ [40, 50] MeV. For the 
remaining inputs My, Ms and F, we consider the same values used before, yielding 
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10 3 L 9 b) 
f 10 s (C 88 -C 90 ) 




ft (GeV) 



Figure 3: The RxT predictions (solid gray band) for the xPT C(p 4 ) low-energy constant 
Lg(fj,) (a) and the 0(p 6 ) combination Css(/j) — Cgo(^) (b) are compared to their large-iVc 
estimates (red dashed) for different values of the renormalization scale fi. The error of the 
large- Nq estimate is given by the naive saturation scale uncertainty from eq. ( |4.9| ). 



the predictions 

L 9 (/i ) = (7.6±0.6)-10- 3 , 
C 88 (/i ) - C 90 (fi ) = (-4.5 ± 0.5) • 10" 5 . (4.12) 

As it can be observed, the influence of using the first or the second approach 
is not crucial at the present level of accuracy. We take the values in ( f4.11| ), which 



include more theoretical constraints, as our final next-to-leading-order estimates for 
the LECs. 

In table |l| we present the different contributions to the LECs within the first and 
second approaches. A graphical comparison of the NLO predictions and the large- A^^ 
estimates has been made in figure |3| for different values of the renormalization scale 
fi. 

It is appropriate to note the appreciable increase of Lg(/io) from the large- Nc 
prediction, for /i = 770 MeV. In fact, the correction S^lo in eq. (|3.6|) gets a contri- 
bution from the An channel which is still comparable to that from the tttt one. This 
subleading contribution to FyGy, fixed through short- distance matching, increases 
the value of L$ by 1 ■ 10~ 3 , a quite sizeable shift. For details see appendix [E] and 



ref. |I|. 



5. Conclusions 

In this article we have completed the analysis of the VFF at NLO in 1/Nc, initiated in 



ref. |TI[ , where the general framework was established. We have considered operators 
with more than one resonance and have studied contributions from intermediate 
channels with resonances. We get a well-behaved VFF at high-energies, which goes 
to zero for q 2 — > oo |26| . 
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10 3 -L 9 (/i ) 


10 5 • (C 88 (/i ) - C 90 (/i )) 


This work 1st 
This work 2nd 
Ref. U 
Ref. [23| 
Ref. [|25j 
Ref. p| at 0(p 4 ) 
Ref. 32 at 0(p 6 ) 
Ref. |53] 


7.9 ±0.4 
7.6 ±0.6 
6.9 ±0.7 
5.93 ± 0.43 
7.04 ± 0.23 
6.54 ±0.15 
5.50 ±0.40 
6.3 ±0.4 


-4.6 ±0.4 
-4.5 ±0.5 

-5.5 ±0.5 



Table 2: Comparison of our result with other determinations, being = 770 MeV. 



Imposing that each individual absorptive cut vanishes at short distances, one 
gets stringent constraints on the structure of the VFF, which led to a prediction 
of the relevant 0(p 4 ) and 0(p 6 ) couplings up to NLO in 1/Nq- The required 
inputs are the resonance masses M v , Ma and M s , and the pion decay constant 
F. As expected for such a well-known observable, the large-iVc prediction provides 
already an excellent estimate and the subleading corrections are relatively small. At 
the reference scale /x — 770 MeV, we obtain 

L 9 (/i ) = (7.9±0.4)-10" 3 , 
C 88 (/i ) - C 90 (// ) = (-4.6 ± 0.4) • 10" 5 . (5.1) 

As the matching of R%T with x?T is complete up to NLO in 1/Nc, we fully control 
the running of the LECs up to that order and, e.g., we are able to predict Lg(fi) for 
any desired value of \i. 

This result is in reasonable agreement with previous calculations 0, E31 BR 



, see table EJ, and shows once more the efficacy of R%T to describe low-energy 
QCD matrix elements, specially if they are dominated by resonances. It is important 
to remark not only that the amplitude is dominated by tree-level exchanges but also 
the fact that the one-loop corrections are not large. 

Our determination of L 9 (/xq) has a larger central value than the result obtained 



from an 0{p ) chiral fit to the VFF |23| at low energies, and it is closer to the 



chiral fit determination at 0(p A ) M. On the other hand, the ALEPH r-data analysis 
performed in f2"5| , which is also of 0(p 6 ) but takes higher-energy data into account, 
yields a value of the order of 7 ■ 10 -3 , much closer to our estimate. 

In future works, we plan to study the pion scalar form-factor and the LECs 
L 4 (/x) and L 5 (fi), where the situation is much less clear since, in that case, one has 
contributions from broad resonance states like the /o(600). 
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A. Dispersion relations and loop contribution 

One may use a once-subtracted dispersion relation, derived from the identity 

m = ±i di _M_ (A.i) 

S 27TZ J t (t - s) ' K ' 

where the integration is performed in the usual complex circuit The form-factor 
in the integrand can be written as 

™ = C "> „ (A.2) 

where D(t) is an analytical function except for the unitarity logarithmic branch cut 
and the single pole of ^p- at t — 0. One gets then 

1 = 1 + 1 J^(.) - Hg^tjgj + RggCj^ f (A . 3) 

where the - term on the r.h.s. is given by the integration ^ fa n v tt^t, with 

s ° J ° 2-ki Jy=2ir~ t (t—s)' 

t = ee ld , around t = of the function m j + O(t ), and the different contributions 
of each two-meson absorptive cut are given by the dispersive integral, 



F U (s)\ mim2 = lim 



s f M v-*lnLF{t)\ mi>m2 , s f°° , lmF(t)\ mi>rn2 



dt / n " n : " 2 + - / dt 

TT JO t(t - S) TT J M 2 +t t (t 



28 lim ((M£-t) 2lmJ7(t)U > 



7TC t-+M£ [ t (t 



(A.4) 



Notice that if the threshold of the channel is above the resonance mass My, then 
this expression gets simplified into the form 

•U/ n \l . i- S f°° j- Ini^wlmi.ma 



^( a )U >ma = hm- / dt ^T , (A.5) 

T J(M 1+ M 2 )2 t[t - S) 
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with Mi (M 2 ) the mass of the mi (m 2 ) meson. 

If we choose the on-shell mass scheme, without double poles in the perturbative 
expansion, we have then 

Ht) = 1+ £ F U (t)\ mi , m2 - SRe ^ v) , (A.6) 

■mi,ni2 V 

where KeD'(M v ) can be identified with — y ^ y - for a convenient renormalization 
scheme of this combination of vector couplings [T7|, |18|, [19j (see appendix [E] for further 
details). 



B. The spectral functions Im J-"(s)| mi) 



m 2 



In this appendix we show the explicit form of the the spectral functions of the different 
two-particle absorptive cuts. First we present the functions obtained directly from 
the Feynman diagrams before imposing any short-distance constraint, i.e., they are 
badly behaved at high energies. 

™» - f ^4-^r H ( mj i ° g ( i + ( - i2mj - es) + s3 

+12sM^j+G 2 v (s 3 -6M v (M 2 +2s) Lgfl+^j (2M 2 + S )-2sj\^ 

s 2 G v (F 2 (F v + 2G V ) (M v - s) + 2sF v G v ) s 

64vrF 6 [s - M 2 ) 2 64vrF 2 ' 1 ' ' 

Im ^^ = 32*F*s\s- M*) Y K (4M| + S) ~ 3M| (2M| + ~ M% 

-6M|(M| - M 2 P ) (-M 2 + 2M| + s )l og (l + ^^) + 12sM| + s 3 j 

+Mf, + hsM 2 + 2M V + 2s 2 ) log(l+ 2 P ) - 3M£ (4M 2 + s ) 
+M* - s (24sM 2 + \2M V + s 2 ) - 2S{S s Zuf } ' (R2) 

lmT(s)\ An = ^ ~ 4 ^~ s M ^ 2)2 {^((2/.+(7)M|+45(K +t T)Ml+ S V) (s-M 2 ) 
-F v (s - Mi) ((2k + <r) 2 Mj + 2s (k 2 + 4cr/« + 2a 2 ) M\ + sV) 
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(kM\ + a(M 2 + s)) + F v ((Ml - s) (M 2 A - M v ) (M v + s)a 2 + 2kM 2 a 
(M 2 A -s)(M 2 A -M 2 )a+K 2 M 2 A (3M A -5(M 2 +s)M 2 A +(M 2 +2s)(2M 2 ^ 
M v +(M 2 A -s) (F A (s-M v ) ((3k+<t)M£+((3<7- 6k)M 2 + s{3k+4ct))M 2 a 
+a (s 2 -6M V - 3s M 2 )) + F v ((M 2 A - s) (M\+AsM 2 A - 6M V + s 2 
+3 (Ml - s) M 2 ) a 2 + QkMI (M 2 a - s) (M 2 A - 2M 2 + s) a 

+k 2 M 2 a (7M A - 8 (3M V + s)M\ + 12M V + s 2 + 24sM v ))) 
V2c d Xf A f , / , s-Ml 



je log(l + (F V (2kMI + (« - a)M 2 



32FHs (s - My) 
(s - Mi) + (k + a)M\ (s - Mi)) + F A (M 2 A - M 2 ) (M 2 - s)) M 2 

+ (Mi - s) (F v (3a (s - Mi) (Mi - 2M| + s) 

+k (4sM 2 A -5M A + 12Ml + S 2 ))+3F A (M 2 A -2M 2 + s)(M 2 -s))^ , (B.3) 

where we have used the combination of couplings k and a, 

k = -2\ V 2 A + Aj A , a = 2A^ A - 2A^ A + Aj A + 2A^ A . (B.4) 

After considering the constraints explained in section 3, 

FyGy = F 2 , 3G 2 v + 2c 2 =F 2 1 

Af v =0, k = 0, (B.5) 

K + CT = ^i a sa_ F A G V {M 2 A -AM 2 ) 

Fv' 1 3V2Mic d F y ' 

the imaginary part of each absorptive cut vanishes at short-distances and the follow- 
ing expressions are found, 

In^( S )U = ^ F ^l_ M , f {™$(F 2 -3G V ) (M* -s)log(l + ^) (2Mf + s) 
+G^M^ log (l + (-6s 3 - 9s 2 M 2 + 6M£ + 9sM*) + 13s 3 

- 6s 2 M 2 - 6sM^j + 6sMl (F 2 - 3G 2 V ) (s — M v ) j , (B.6) 

ImF(s)\ Pn =0, (B.7) 
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V h 32vFHM\(s-Mlf\ A \ SK V J \ V 



i- A ys-iviy) i. \ \ ■ Mj ' 

+4 S M 2 - 7My - 3s 2 ") + 2M 2 ( S 2 M y ( 3 logfl + i^S) - 2 



1 



My 



+Mj ( a - My) log(l + - Ml (5 - My) 



+M V (s*M v (l- 6 log (l + ) + 8MI (M v - s) log (l+ iJjp) 

+6M y log (l + ^r^) +8sMj (s - M v ) -QsM^j + 2sM\ + M* j . (B.8) 



C. Next-to- leading-order corrections ^ NLO (s)| 



NLoV* J /|mi,m 2 



In this appendix we show the explicit form of the NLO corrections generated by the 
considered two-particle absorptive cuts, eqs. ( P-6| ), ( p.7j ) and ( P-8|) , which have been 



calculated by using the dispersive method discussed in appendix Below, we have 
summed up the <5 NLO contribution to J rU (s), as seen in eq. (|3.6|) , being the different 
J r NLO (s)| mim2 well-behaved at high energies: 

^Wl« = u^s-mA - 12Ml (F2 - 3Gl) {s - Ml) /(s ' M|) 

- 6MI (F 2 - 3G V ) (s - M v ) (sf( 8 , M 2 ) + 2 log (^J) - 2 ) 



+ G v My ( - 6 (3s 2 M v - 3sM v - 2M V + 2s 3 ) f(s, M v ) 
+s2 (- 26l0g (Mf) +27 ) +12M K l0g (5)" 1 ) 



(C.l) 



+ 3.M 2 ( 4 log (^) - 5)) + 3sM 2 s (F 2 - 3G V ) (s - My) 
F NLO (s)\ Pn = , (C.2) 

^ - -3^S^{ m ^H w - s) 9 (s « 

-Qs 2 log(l - ||) + log(l - ^) (3, 2 +2M 2 (M 2 - S ) + 7M v -4sM$) 



+ sM v [My -6 (s 2 - My) g (s, M 2 , M y ) + 6M 2 log 1 



M 2 
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+8Mj(M* - ,),(,, Ml, MJ) -8MJ( S - Mj) (log(l- JL) + l g(^|) 
-l)) + M\M' V {m 2 U (s 3 - sM' v ) g (s, Ml Ml) + s 2 (a log (l - JL) 

+2bg(l-g)-7)-«Mjlog(l-^) + 7^) 
+2sM|( S -M^)^( S) Mi,M|) 

+2M|( S - My) (4 S M^ 5 ( a , Mi, M|) + S (log(l - ^) + log(^) - l) 
+ 4M^ log (l - ^) ) ) + Mi log (l - g) - M* log (l - ^) ) 
+ sM^My (m 2 (-2 log (l - ^) " l) + 2^ log(l - g) + } , (C.3) 
where the functions f(s,M 2 ) and g(s, M±, M$) have been introduced for simplicity, 

'(••"*> = §H 1 + 

9( .,«?,^) = i(u 1 (i + ^-g)-i J .(i-g)). 



(C.4) 



D. NLO contributions to Lg(/i) and Cssi/j) — C9o(m) 

In this appendix we give the full expressions of the NLO contributions to Lg(n) and 
C 88 (/i) — Cg (fi), following the notation of eqs. ( P~op and Q£77fl, i.e., £ml,m 2 and £mi,m 2 : 

® ^ SiM* - Mm}) - ") + ( 38 - 181 -(m|)) } ' M 

62 = , (D.2) 

t(2) _ F A G v /2M 10 (M 2 - M 2 ) lotrfl - 

U - l 287r 2 F 4 M 2 M 8 (M 2 _ M |) |2M A LM fi M A ) log^l M , j 

-2M1M 2 (Mi - Ml) (log(l -g)+l)+ Mi My (Mi - Mf) 
(l61o g (l-g)-3) + M|My(M|(-21o g (g) + 161o g (l-g)-ll) 
+ Mi (ll-161og(l-§ ))) + MiMy (MJ (l01og(f )-121og(f ) 
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- 2 .o g ( 1 -f) + n) + M,( 121 o g (f) +2log ( 1 -f)-n)) 

+ M-(Ml(-61o g (g) +2 .o g (l-g) +5 ) 

-^( 81 -(i)- 61 -(8) +2l0g ( l -m) +s ))l (D3) 

«- (» - m®) + *(»*(3) - ") } 

S£ = 5 1 - 6MA°M2 (Mi - M|) 2 

A " 3847r 2 F 2 M|M^° (M| - M|) (M| - M v ) \ A VK A S> 

+ 8M- {Ml (log(g) + l) -M^ -M- (2MiM| ( 3 tog (§) +6 ]og(g) 

, / M 2 x \ ,/ , /M^\ , /M 2 x , / M^\ \ 

+AfJ (-61og(M) +6 l og (l _ ) _ 13^ _ SMlM^^MlMj 

(^(f)---0)-f-<-^(i)— 



+ M1 (61o g (^) - 9 )) + 2A« (m^M| ( 8 log(g) - 18 log(^) 

(91«g(^) + 27 log (l - M) _ 35 )) - M\Ml (M\ - Ml) 
(Ml (3 log(g) -96 log(l - g) + 4 7 ) + Mi (96 log(l - g) -4 7 )) 
+3Af» (Afl - M| 2 (l8 log (l - _ 1] _ 6 Af- (Mj - M| 2 log(l - ^) 



(D.6) 



E. Description in terms of Feynman diagrams 

The subleading corrections can be calculated by means of dispersive relations. Once 
the NLO absorptive parts of J-r x t (s) are known, one can reconstruct the full form 
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factor up to appropriate subtraction terms. Alternatively, we can compute and sep- 
arate the tree-level and one-loop amplitudes in the form 



_ FyGy S 2Zvg ^ — \ 

Fr-xT (a) = 1 + f2 m2 _ ~ + -pT s + 2^ ^(^"w (E- 1 ) 

^ mi,rrt2 

where the one-loop diagrams J r (s)| mi)m2 can be rewritten by means of a once- 
subtracted dispersion relation in the form 

E_, , | ^ > / \ i 2(5^2 r ^ ? ^ 

■^( s )|mi,ma = ^ l S )Ui,TO2 + ^2 S + "0^2 _ ~ + d -2/^2 _ g \2 ■ 

rni : rn '2 m i ,mj ^ \ v / 

(E.2) 

The finite part of the loops is contained in the once-subtracted dispersive func- 
tions J~ u {s)\ mum2 , fully determined by the imaginary part of ImJ r (s)| mi)ma through 
eq. ( |A.4|) . The real parameters £-2,0,2 contain the ultraviolet divergences of the loops, 
being <5 and 5-2 the real part of the pole residues. The local R^T coupling Lg renor- 
malizes 62 , the combination FyGy cancels the divergences in 5q and a convenient 
shift of the mass, jl4 B) 2 = My + SMy removes the divergent part of 8-2- Indeed, 
we will work in the on-shell scheme and the counterterm 5M V will be chosen to 
completely kill 6-2- 

In order to finish the short- distance matching we just need to take into account 
that the once-subtracted loop contribution behaves at short distances like 

mi,rn,2 

with 5o a constant number (denoted before in the text as S NLO ). This leads to the 
VFF high-energy constraints 

—=2 h o = 1 + o , 

F 

L 9 + 5 2 = 0. (E.4) 

Hence, the VFF finally takes the well-behaved structure (|3.5|) employed in the 
article, 

s 



, r ' / , ~ ( S ) I"' I '"'J 

My — S 

v mi,rri2 



My 



S/f2 +^ NLO (^)- (E.5) 

My — S 



Notice that no real double pole term 5_ 2 remains in our perturbative NLO expression 
as we have chosen the on-shell mass scheme. 
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